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We study chaotic dynamics in preheating after inflation in which an inflaton (f> is coupled to 
another scalar field % through an interaction (l/2)g 2 (/) 2 x 2 ■ We first estimate the size of the quasi- 
homogeneous field x at the beginning of reheating for large-field inflaton potentials V(<fi) = Vo4> n 
by evaluating the amplitude of the \ fluctuations on scales larger than the Hubble radius at the end 
of inflation. Parametric excitations of the field x during preheating can give rise to chaos between 
two dynamical scalar fields. For the quartic potential (n = 4, Vo = A/4) chaos actually occurs 
for g 2 /X < O(10) in a linear regime before which the backreaction of created particles becomes 
important. This analysis is supported by several different criteria for the existence of chaos. For the 
quadratic potential (n — 2) the signature of chaos is not found by the time at which the backreaction 
begins to work, similar to the case of the quartic potential with g 2 /A 3> 1. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Reheating after inflation is an extremely important 
stage to generate elementary particles present in cur- 
rent universe. In the original version of the reheating 
scenario which is now called old reheating, the decay of 
an inflaton field is characterized by a perturbativc Born 
process 0. However this process is not efficient for the 
success of the GUT-scale baryogenesis scenario. Later it 
was found that the existence of a nonperturbative stage- 
dubbed preheating- can lead to an explosive particle pro- 
duction prior to the Born decay 0- B Wl 

During preheating scalar particles x coupled to the in- 
flaton (f> are efficiently generated by parametric resonance 
through an interaction (l/2)g 2 cf> 2 x 2 - The existence of 
the preheating stage provides several interesting possibil- 
ities such as the GUT-scale baryogenesis ["} , nonthermal 
phase transition |~j , the enhancement of metric perturba- 
tions and the formation of primordial black holes 
In the chaotic inflationary scenario characterized by the 
potential V(cj)) = Vo<j) n , the field perturbations <5% obey 
the Mathieu equation (for n = 2) or the Lame equation 
(for n — 4), which determines the structure of resonance 
at the linear regime. When the backreaction of created 
particles begins to violate the coherent oscillation of (f>, 
this tends to work to suppress exponential growth of the 
field fluctuations. The system enters a fully nonlinear 
stage after which the mode-mode coupling (rescattering) 
between perturbations is crucially important [tj, ITol | . 

Typically the contribution of the dynamical back- 
ground field x is neglected in standard analysis of parti- 
cle creations in preheating. This may be justified for a 
quadratic inflaton potential, since large-scale x modes are 
exponentially suppressed during inflation for the coupling 
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g required for preheating [Tlj . However the situation is 
different for a quartic inflaton potential with the coupling 
g of order g 2 /X = 0(1) fl"l Il3|. In this case the quasi- 
homogeneous field x can pl a Y an important role for the 
dynamics of preheating. In fact Podolsky and Starobin- 
sky pointed out that chaos may occur for the self- 
coupling potential V(4>) = (1/4) A0 4 when the coupling 
g 2 /X is not too much larger than of order unity. Since 
it is not obvious whether chaos actually occurs or not 
in this model only by analytic estimations, we shall per- 
form detailed numerical investigation with/ without the 
backreaction effect of created particles. We shall esti- 
mate the size of the field x a t t ne beginning of reheating 
by evaluating the amplitude of the x fluctuations for the 
modes larger than the Hubble radius. To judge the ex- 
istence of chaos in preheating we will adopt several dif- 
ferent methods- such as the Toda-Brumer test 0, ^1 , 
Lyapunov exponents |l7j | and a fractal map [i"~f . 

It was already found that chaos appears for hybrid- 
tvpe inflation models [H HO, |2~1 |~~| (see also Refs. |H 
|2J, |25j). Hybrid inflation is a rather special model in 
a sense that the symmetry breaking field automatically 
grows by tachyonic instability even if it is suppressed dur- 
ing inflation. The necessary condition for chaos is that 
there exist at least two dynamical fields and neither of 
them is too much smaller than another field. The hy- 
brid model satisfies this condition, since two fields can 
have frequencies which are the same order after symme- 
try breaking. The presence of mixing terms between two 
fields leads to a new instability of perturbations in addi- 
tion to tachyonic/resonance instabilities |2jJ. This new 
type of instability is clearly associated with the presence 
of chaos. 

In this work we shall investigate the existence of chaos 
for large-field potentials V(4>) = Vo<f> n with an interac- 
tion (\/2)g 2 (j) 2 x 2 ■ We estimate the variance of large-scale 
modes in x a t t ne end of inflation, which is relevant to 
the initial condition of the quasi-homogeneous field x f° r 
preheating. The field x IS amplified by parametric reso- 
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nance, which can give rise to chaos after X grows to satisfy 
the Toda-Brumcr condition. We shall numerically solve 
background equations together with perturbed equations 
for both quartic (n = 4) and quadratic (n = 2) poten- 
tials. Our main interest is to find the signature of chaos 
and the parameter range of the coupling g in which chaos 
can be seen before the backreaction effect of created par- 
ticles becomes important. Since chaos can alter the stan- 
dard picture of preheating by parametric resonance, it is 
of interest to clarify the situation in which chaos appears. 

Recent observations suggest that the quartic poten- 
tial V(4>) = (l/4)A</> 4 is under an observational pressure, 
while the quadratic potential V(<p) — (l/2)m 2 2 is al- 
lowed j2|| . This depends on the number of e-folds TV 
before the end of inflation at which observable pertur- 
bations are generated. In the case of quartic potential 
this corresponds to N ~ 64 by assuming instant tran- 
sitions between several cosmological epochs [27|. The 
likelihood analysis including WMAP and SDSS datasets 
shows that the quartic potential is marginally allowed 
by using N = 64 p3 |. Therefore it is premature to 
rule out this model completely from current observations. 
The quartic potential corresponds to a system in which 
the background equations can be reduced to the ones in 
Minkowski spacetime by introducing conformal variables. 
This has an advantage for the investigation of chaotic dy- 
namics during preheating. As we see later, the quartic 
potential exhibits a stronger chaos compared to the one 
for the quadratic potential. 



II. THE FIELD VARIANCE FOR LONG 
WAVELENGTH MODES 



We define the number of e-folds before the end of in- 
flation, as 



N = \n(a f /a(t)) , 



(5) 



where a/ is the scale factor at the end of inflation. Em- 
ploying the slow-roll approximation, \(j>\ <C |3Zf</>| and 
4> 2 <C V{4>), we easily find that d^/diV = n/(n 2 (f>). Here 
we neglect the contribution coming from the x-dependent 
terms. Integrating this equation, we obtain 



2n 



-N. 



(6) 



The field value at the end of inflation ((f) f) is determined 
by setting the slow-roll parameter, e s = (1/2k 2 )(V<p/V) 2 , 
is unity. This gives 4>f/m p = n/(4y / 7r), thereby leading 
to 



= TA N+ i)< 



(7) 



Let us consider a perturbation Sx in the field \. Then 
the momentum-space first-order perturbed equation is 
given by 



s Xk + 3H6 X k + H> + srr &xk = o , 



(8) 



where k is a comoving wavenumber. Note that we ne- 
glected the backreaction of gravitational perturbations. 
When the effective mass of the field X 1S larger than of 
order the Hubble rate, i.e., g 2 (f> 2 > (3H/2) 2 , the evolu- 
tion of super-Hubble perturbations is characterized by 
underdamped oscillations with the dependence 



In this section we shall estimate the field variance for 
long wavelength modes of the field x coupled to the infla- 
ton <f> through an interaction (l/2)g 2 cf) 2 x 2 - The effective 
potential in our system is 



v(<t>,x) = Vo4> n 



1 
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We are mainly interested in two inflaton potentials: (i) 
the quadratic one (n = 2) and (ii) the quartic one 
(?i = 4). In this work we do not implement nonmin- 
imal couplings [2^| between the field x an d the scalar 
curvature R. 

In a flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) background with a scale factor a, the back- 
ground equations for the system Q are 



(f> + ZH4> + = , 
X + 3H X + V X = 0, 



H 



where Va, 
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(a, 


) -t[ 



\<t> 2 + \x 2 + v(4,x) 



(2) 
(3) 

(4) 



V x = g 4> x an d k 



8Tr/m 2 with m p being the Planck mass. 



5 X k oc a 3/2 



(9) 



Hence large-scale perturbations are exponentially sup- 
pressed during inflation. Meanwhile when g 2 (p 2 < 
(3H/2) 2 super-Hubble perturbations evolve as 



Sxk oc exp 




(10) 



This shows that in the massless limit (g 2 <fi 2 <IC H 2 ) the 
amplitude of Sxk decreases very slowly. In what follows 
we shall consider the quadratic and quartic models sep- 
arately. 



A. Quadratic model 



For the inflation potential V{4>) = jin f, the coupling 
g is required to be greater than of order 1CP 4 in order 
for preheating to occur [T^j. In this case the resonance 
parameter, q = q 2 4> 2 /(Am 2 ), is much larger than 1 at the 
beginning of preheating |l0j . Since H ~ m at the end of 
inflation, the existence of the preheating stage demands 
the condition g 2 (f> 2 ^> H 2 during inflation. Hence the 
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evolution of large-scale x fluctuations is characterized by 
Eq. ©. 

Let us consider the modes which are outside the Hub- 
ble radius at the end of inflation (0 < k < kf — afHf). 
Since these modes are effectively massive with slowly 
changing mass, they can be treated as an adiabatic state 
characterized by 5\k = l/(a 3 ^ 2 \/2uJk) where lo\ ~ g 2 4> 2 - 
Hence the amplitude |<5xfc| 2 at the end of inflation is es- 
timated as 



\6Xk(t f )f 



1 



2a a f gfa 



(11) 



Then we can obtain the variance of the fluctuation Sxk 
for k < kf-. 

(Sxl(tf))k<k f = ^ f ' k 2 \6 Xk (t f )\ 2 dk 



(12) 



where we used </>/ = m J) /(2y / 7r) in addition to the slow- 
roll approximation H 2 ~ 47rm 2 2 /3m 2 . 

For the quadratic potential (n = 2) the inflaton mass is 
constrained to be m ~ 10~ 6 m p from the COBE normal- 
ization j^J ■ Then Eq. I|12|l shows that the variance is the 

function of g only. In Fig-.^we plot J (#X/b(*/))fc<fc/ as a 
function of g. It can be regarded as a quasi- homogeneous 
mode at the beginning of preheating. This at least mea- 
sures the minimum amplitude of the homogeneous \ field. 




B. Quartic model 

For the quartic potential V(4>) = (1/4) A<fi 4 , it is known 
that preheating occurs even when the coupling g is in the 



range g 2 2 < H 2 when 9 2 / x = for exam- 

ple, the background dynamics transits from the "mass- 
less regime" [g 2 <i> 2 < (3H/2) 2 ) to the "massive regime" 
[g 2 4> 2 > (3ff/2) 2 ] during inflation EHHHI. The crit- 
ical number of e-folds, N c , in which the the evolution of 
the perturbation 8\k transits from Eq. i|l(J|) to Eq. @ is 
determined by the condition 9H 2 /4 = g 2 (j> 2 , which gives 

M 



Nr. 



= In 



In 



3 A \a c J \k c 

We note that we used the slow-roll condition 

2n\(j) 4 ,„ N o 



H 2 



3m 2 , 



-mz 



(13) 



(14) 



The modes which are inside the Hubble radius at tran- 
sition time t ~ t c (but larger than the Hubble radius at 
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FIG. 1: The amplitude of the variance {<5Xk(~tf))H<k f a ^ the 
end of inflation in terms of the function of the coupling g for 
the quadratic potential (n = 2). The variance gets smaller 
for larger g. 



the end of inflation) evolves as effective massive fields for 
t > t c . Then for k c < k < kf the amplitude |<5xfc| 2 at 
the end of inflation is given by Eq. (|llf> . which gives the 
variance 



(Sxl(tf))k c <k<k f 




k 2 \S X k(tf)\ 2 dk 



)m 2 



(15) 



where we used the relation k c = kf exp(— 2g 2 /3A) coming 
from Eq. JT^J. 

For k < k c the modes exit the Hubble radius before 
the transition time t = t c . At the Hubble radius crossing 
(t = tk) the amplitude of the perturbation Sxk is given 
by 



\5 X k{tkt 



H 2 {t k ) 
2k 3 



(16) 



which comes from the quantization of a standard massless 
scalar field [3(3]. Since 8\k evolve as Eq. I|1(J|I for t k < 
t < t c , we find that the amplitude of the perturbations 
at t = t c is given by |33l ] 



\6Xk(t f )f 



H 2 (t k ) 
2fc 3 ' 



-3F(N k 



(17) 
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where 



F(N k ) = N k -N e - VN k (N k - N c ) 
m~ k + ^N k -N c 



-AUn 



(18) 



Here JVfc is the number of e- folds at t = t k . The pertur- 
bations evolve as Eq. ljTU|) for t c < t < tf. Using Eq. ijT3|) 
we obtain the following amplitude at the end of inflation: 



_ H 2 (t k ) n - 3 F(N k )-2g 2 /\ 



\sxi(t f )\ 



2fc 3 



(19) 



Then the variance of perturbations 8\k for the modes 
k < k c is given by 



(SxUt f )h< kc 



1 

2^ 



ISxkitfWkMilnk) 



k, 

1 



\ rl ,2 r Ni 

= j dAWfe" 3 ™" 2 ^ ,(20) 

where we used Eq. <|14|) and d(lnfc) ~ d(lna) = —dN k . 
Note that hi is the minimum wavenumber relevant for 
the maximum scale of cosmological perturbations. In 
order to obtain {8x\{tf))k<k c , it is necessary to solve the 
following differential equation: 



dN, 



-<*Xfc(«/))fc<fc. 



Xm 2 

6-7T 3 



N 2 k e 



-3F{N k )-2g 2 /X 



(21) 



which should be integrated from N c = 2g 2 /3X to Ni. 
Note that Ni roughly corresponds to the total number of 
e-folds during inflation. At least we require the condition 
Ni > 60. We shall choose Ni = 60, 100, 1000 in order to 
see the sensitivity for the change of this number. 

Finally the variance for the modes k < kf is given by 
the following sum: 

(Sxl(tf))k<k f = (Sxl(tf))k c <k<k f + (Sxl(tf))k<k c -(22) 

When g 2 /\ < 0(1) one has N c < 0(1). This shows that 
most of the contribution to the variance (SXk(tf))k<kt 
comes from the modes k < k c . Meanwhile when g 2 /X 3> 
1, i.e., JV C > 1, the modes k c < k < kf dominate the 

total variance. In Fig. |21 we plot J (Sx k (t /)) k< k f as a 

function of g 2 /\ for N t = 60, 100, 1000 and A = 10~ 13 . 
We find that the variance does not depend on the values 
Ni for g 2 /\ > 3, which reflects the fact that Eq. lfT3|) is 
independent of N. The difference appears for g 2 /X ^ 3, 
because of the fact that (5x k (tf))k<k c is dependent on 
the values N. We shall use the values (Sx k (tf))k<k f ob- 
tained for Ni = 60 as a minimum initial condition of the 
quasi-homogeneous x held at the beginning of preheat- 
ing. 

We note that Podolsky and Starobinsky [lj| estimated 
the size of the quasi- homogeneous x ne ld f° r ff 2 M = 
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FIG. 2: The amplitude of the variance (8x%{tf))]/< k at the 
end of inflation in terms of the function of the coupling g for 
the quartic potential (n = 4). The variance is dominated by 
the modes k c < k < kf for g 2 /\ 3> 1, whereas the dominant 
contribution comes from the modes k < k c for g 2 /\ < 0(1). 



by using a Fokker Planck equation |3J| . This equation is 
valid when the mass of the field x is smaller than of order 
the Hubble rate H^,|3^,|33, which corresponds to g 2 /X ^ 
0(1). We checked that our estimation of the variance 
shows good agreement with the one based on the Fokker 
Planck approach for g 2 /X < 0(1)- In the parameter 
regime g 2 /\^> 1 the Fokker Planck equation is no longer 
valid. Hence we need to use our estimation given above 
in order to know the size of the quasi-homogeneous x 
field at the end of inflation. 



III. BASIC PROPERTIES OF PREHEATING 
AND CHAOS 

A. Preheating and the role of the 
quasi-homogeneous field \ 

In the presence of the coupling (l/2)g 2 cj) 2 x 2 , the co- 
herent oscillation of the inflaton leads to the excita- 
tion of the field x during preheating through the reso- 
nance term g 2 (j) 2 x in Eq. For the quadratic potential 
(n = 2) parametric resonance is efficient when the con- 
dition, q = g 2 cf> 2 /(4m 2 ) 3> 1, is satisfied. To be more 
precise the field x grows for g > 3.0 x 10~ 4 by over- 
coming the friction due to cosmic expansion [Tol | . In this 
model the growth of x ends when the system enters a 
narrow resonance regime [q < 1) or the backreaction ef- 
fect of created particles breaks the coherent oscillation 
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of the field </>. For the quartic potential (n = 4) reso- 
nance bands exist for the parameter space characterized 
by n(2 n - 1) < g 2 /X < n(2n + 1), where n = 1,2, 3, ... 
|3lLl32l |. The center of the band, g 2 /X = 2n 2 , corresponds 
to the largest Floquet index (/i max — 0.28). 

If the field x is strongly suppressed during inflation, 
this does not contribute to the background dynamics 
even if it is amplified during reheating 1 . As we see in 
Figs. H an d El it is expected that x does not dynami- 
cally become important during preheating for larger g. 
However, for the quartic potential, parametric resonance 
takes place even for g 2 /X = 0(1), in which case the quasi- 
homogeneous x ne ld does not suffer from strong sup- 
pression during inflation. Then if the quasi-homogeneous 
field x is amplified during preheating, this can contribute 
to the background dynamics. 

In our model the equations for field perturbations in 
Fourier space may be written as 0] 



b k + 3H8cf> k + 
-2g 2 (j>x5xk , 



+ n{n - l)V^ n - 2 + g 2 X 2 



5 Xk + 3H5 Xk + [ — + g 2 4> 2 ) S Xk 



-25' 



0(pk 

(23) 
(24) 



where k is a comoving wavenumber. This corresponds 
to the equations in which the contributions from metric 
perturbations are dropped (see Ref. Q). Since metric 
perturbations are enhanced only when field fluctuations 
grow sufficiently, it is a good approximation to neglect 
them except for a nonlinear stage of preheating. 

The terms on the r.h.s. of Eqs. (|23|l and (|24|1 are not 
usually taken into account in standard analysis of pre- 
heating 0, 0> E3 > since the field X was supposed to be 
dynamically unimportant. However this can give consid- 
erable contributions to the dynamics of field perturba- 
tions provided that the quasi-homogeneous field x is not 
strongly suppressed during inflation. In fact these terms 
lead to a mixing between the perturbations of two fields, 
whose behavior is absent in standard analysis of preheat- 
ing unless rescattering effects are taken into account at a 
nonlinear stage. 

We start integrating background and perturbation 
equations from the end of inflation. We adopt the Bunch- 
Davies vacuum state for the initial condition of perturba- 
tions for the modes inside the Hubble radius. The total 
variances of the fields ip = <p, x integrated in terms of k 
are 



(^ 2 ) = ^2 / k 2 \5tp k \ 2 dk. 



(25) 



1 The amplification of \ is limited by the backreaction effect of 
created particles. 



We implement the variances (5<p 2 ) and (5x 2 ) for both 
background and perturbation equations as a Hartree ap- 
proximation We note that this is for estimating 
the time at which backreaction effects become impor- 
tant. After the system enters a fully nonlinear stage, one 
can not trust the analysis using the Hartree approxima- 
tion. Our interest is to find a signature of chaos before 
the backreaction sets in. 



B. The condition for chaos 

In this subsection we review several conditions for the 
existence of chaos and apply them to our effective po- 
tential QJ. Let us consider the first-order differential 
equations 



X{ F-i (xj 
and their linearized equations, 
- . BR 



(26) 



(27) 



where Sxi is the perturbation vector connecting two 
nearby trajectories and dFi/dxj is the Jacobian matrix 
of Fi(xj). 

The scalar-field equations J3J) and are expressed by 
the form l|2tj|) by setting x\ — 0, xi = 0, £3 = x an d £4 = 
X- Then one can evaluate the Jacobian matrix dFi/dxj 
and its eigenvalues [i for a general system characterized 
by an effective potential V — V{<j),x)- Note that we 
do not account for the linearized equation for H, since 
metric perturbations are neglected in our analysis. The 
eigenvalues of the matrix dFi / dxj are given by 



"=2 



-3H ± yjm 2 + 4 7 



(28) 



where 



+ v xx ? 



(29) 



The necessary condition for the existence of chaos is 
that one of the eigenvalues is at least positive. In an 
expanding background (H > 0) this corresponds to the 
condition 7 > from Eq. Ij28(l . Since we are now consid- 
ering a situation in which both effective masses of <j> and 
X arc positive (V^ > 0, V xx > 0), 7 can take a positive 
value when 



VL < . 



"0 X ^ u • (30) 

This is so-called the Toda-Brumer test 0, that is 
used to judge the existence of chaos. For our effective 
potential Q this translates into the condition 



x 2 > 



n(n — l)Vb 
3? ' 



(31) 
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When n = 2 and Vq = m 2 /2, this corresponds to 

m 



X > 



whereas for n — 4 and Vq = A/4 we get 



(32) 




(33) 



For the quadratic potential the initial value of \ is 
10~ 4 \/ff times smaller than the value which leads to 
chaotic instability, see Eqs. (|12|l and (|32|l . Therefore the 
system is expected to enter a chaotic phase after the field 
is amplified more than l0 4 /yfg times. For the quartic po- 
tential the condition for chaos is not so severe compared 
to the quadratic potential, since the term on the r.h.s. of 
Eq. l(3*3*|l decreases with time. 

It is worth commenting on the difference about the in- 
stabilities of chaos and parametric resonance. Although 
the field x exhibits an exponential growth by resonance, 
this is different from the chaotic instability in which the 
evolution of the system is very sensitive to slight change 
of initial conditions. In fact none of the eigenvalues of the 
Jacobi matrix is positive in the regime where the condi- 
tion 1)3 1J) is not satisfied. This means that chaos is absent 
in the region \ 2 < n(n — 1)Vq/ (3g 2 )cf> n ~ 2 even if the field 
shows an exponential growth by parametric resonance. 

Since the Toda-Brumer test is not a sufficient condi- 
tion for the existence of chaos, we shall use other criteria 
as well such as Lyapunov exponents ^3] ■ The Lyapunov 
exponents measure the logarithm of of the expansion of 
a small volume in a ^-dimensional phase space. In this 
case the system possesses N Lyapunov exponents. Chaos 
is accompanied by an increase of the size of the N- volume 
at least in one direction and a maximal Lyapunov expo- 
nent h characterizes the signature of chaos. When h ap- 
proaches a positive constant asymptotically, this shows 
the existence of chaos since the initial displacement of 
the iV-volume grows exponentially. If h approaches to 
asymptotically, this means the absence of chaos since the 
orbits are periodic or quasi-periodic. Strictly speaking 
Lyapunov exponents are defined in the limit t — > oo. Nev- 
ertheless one can check the existence of chaos by inves- 
tigating the behavior of the system for sufficiently large 
values of t. In fact, as we see later, the maximal Lya- 
punov exponent begins to grow toward a constant value 
when chaos appears. Although the backreaction effect of 
created particles can alter the background dynamics, it 
is possible to see the signature of chaos before the back- 
reaction sets in. 

In addition to the above two criteria, there exists an- 
other criterion for the existence of chaos- which is so 
called a fractal map 0] . This strategy is useful because 
of gauge independence which comes from using a topolog- 
ical character. In the next section we shall also use this 
criterion to confirm the presence of chaos in addition to 
other methods. 
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$ 

FIG. 3: A phase-space trajectory in the ((/>, x) plane for 
g 2 /\ = 2. The dotted curve corresponds to the boundary 
given in Eq. <l3"?t . 



IV. CHAOTIC DYNAMICS FOR THE QUARTIC 
POTENTIAL 

For the quartic potential (n = 4) the system is ef- 
fectively reduced to a Hamiltonian system in Minkowski 
spacctime by introducing conformal variables: (f> = atp 
X = ax and rj = J a~ l dt. Using the fact that a oc r\ 
and a! /a, a" /a — ► during reheating in this model, the 
background equations © and J3J can be written as 



+ A</> 3 + ff 2 x 2 </> = 0, 



x" + g 2 4> 2 x 



8tt 
3m 2 

8tt 
3m 2 



1 - 



12 



0, 



9 72 ~2 

y0X 



E = const , 



(34) 
(35) 



(36) 



where a prime denotes the derivative with respect to con- 
formal time 77. The trajectories of the fields are bounded 
by 



A 



\~^X 2 <E. 



(37) 



This boundary is plotted as a dotted curve in Fig. |3J 

From the Toda-Brumer test l|33ll . we can expect that 
chaos occurs when x becomes comparable to (j> for g 2 /A = 
O(l). When 1 < oVA < 3, corresponding to the first res- 
onance band [3ll l32l | , the variance of the field x is larger 
than Xf ~ 10~ 8 TO p right after the end of inflation from 
Fig. [3 We require the parametric excitation of x to give 
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x 



FIG. 4: Evolution of X , {$X 2 ) 1/2 and {S4> 2 ) 1/2 (normalized 
by m v ) as a function of x = \/~\4>iri for ,g 2 /A = 2 when the 
backreaction effect of created particles is neglected. The hor- 
izontal line shows the border of the Toda-Brumer test. The 
perturbations in <j> and \ exhibit instabilities associated with 
chaos once the field x is sufficiently amplified. 
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FIG. 5: Evolution of the maximal Lyapunov exponent h for 
g 2 /A = 2, 8, 18. We do not implement the backreaction effect 
of created particles. We find that the exponent h approaches 
a constant value in these cases. 



rise to chaos, since x is much smaller than <fi at the begin- 
ning of reheating. Therefore the coupling g needs to lie 
inside of the resonance band for the existence of chaos. In 
Fig. |3 we show an example of the background trajectory 
for g 2 /X = 2. We find that two fields evolve chaotically 
in the phase-space of the (</>, x) plane by choosing several 
different initial conditions. 

When g 2 /X = 2 the Toda-Brumer test gives the con- 
dition x/rrip > 0.5 for the existence of chaos, which cor- 
responds to the time x = s/\<f>ir) > 45 (see Fig. 0}. We 
find that the fluctuations (Sx 2 ) and (5cf> 2 } exhibit rapid 
increase with a similar growth rate after the field x satis- 
fies the Toda-Brumer test. It is expected that this is asso- 
ciated with the presence of chaos rather than parametric 
excitation of the x fluctuation. The quasi-homogeneous 
field x is amplified by parametric resonance for x < 45 
but stops growing after that. This comes from the fact 
that the resonance does not occur once the homogeneous 
oscillation of <j> is broken by the growth of x- The g cj) 2 
term on the l.h.s. of Eq. (|24|) also becomes ineffective at 
this stage, but the presence of the mixing term on the 
r.h.s. leads to a new type of instability associated with 
chaos. The rapid growth of (Sep 2 ) seen in Fig.^Jalso comes 
from the mixing term on the r.h.s. of Eq. H23[) rather than 
from the parametric excitation of sub-Hubble modes with 
3/2 < k 2 /{\(j) 2 ) < s/3 0,113. 

In order to check the existence of chaos, we plot the 
evolution of the maximal Lyapunov exponent h for sev- 
eral different values of g in Fig. |SJ The exponent de- 
creases at the initial stage even though x is amplified 



by parametric resonance. However h begins to grow af- 
ter x satisfies the Toda-Brumer test (|33f) . The system 
eventually approaches a phase with a positive constant 
h, which shows the existence of chaos. The growth of the 
maximal Lyapunov exponent is regarded as a signature 
of chaos, since this behavior can not be seen in the ab- 
sence of chaos. We checked that h continues to decrease 
and converges toward in power if the mixing terms do 
not exist on the r.h.s. of Eqs. (J2HJ and pi)). 

In Fig. El we show a fractal map for g 2 /X = 2 with 
slight change of initial conditions in terms of xi and xi- 
We set exit pockets when the field x becomes larger than 
|x| = 1.2. When an orbit reaches a pocket we assign 
colors to many initial conditions as in the following way; 
white if an orbit falls down to an upper pocket (x > 1-2), 
black if it falls down to a lower pocket (x < —1.2). Figure 
is the result of the above manipulation, which shows 
that the map of initial conditions is fractal. This means 
that orbits are sensitive to initial data, thereby showing 
the existence of chaos. 

The above discussion neglects the backreaction effect 
of created particles. If we account for it as a Hartree 
approximation, Eqs. 10), ©, Q, HSJ and (JUJ) are mod- 
ified by replacing the 4> 2 , x 2 and (ft 3 terms for (j> 2 + (S(j) 2 ), 
X 2 +(Sx 2 ) and 3 +30(<^ 2 ), respectively [lij. By Eq. d| 
the backreaction becomes important when 

VW)Zj^<t>. (38) 



7.4 x 1(T 7 8.0 x 1(T 7 8.6 x 1(T 7 




7.41 x 10~ 7 7.43 x 10~ 7 7.47 x 10~ 7 

X, 



FIG. 6: A fractal map for g 2 /\ = 2. We show the map of 
initial conditions \i an d XI with exit pockets characterized 
by Ixl > 1-2. The white color corresponds to orbits which 
give x > 1-2, whereas the black one to orbits which give 
X < —1.2. The upper panel corresponds to the change of 
initial conditions by 0.1%. The lower panel is an extended 
figure, in which initial conditions change by 0.005%. These 
figures exhibit fractal structures, which result from sensitivity 
to initial conditions. 



This is similar to the necessary condition for chaos for the 
background field x, see Eq. When g 2 /X = 0(1), 

the Toda-Brumer test is satisfied before the condition 
(J2HJ) is fulfilled. As seen in Fig. U| both x 2 and (Sx 2 ) 
have similar amplitudes initially, but the growth of sub- 
Hubblc fluctuations occurs later than that of x 2 ■ At the 
time when the Toda-Brumer test l|33[) is satisfied (x ~ 
45), (Sx 2 ) is much smaller than x 2 for g 2 /X = 2. The 
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FIG. 7: Evolution of x, (Sx 2 ) 1 ' 2 and (54> 2 } 1/2 (normalized by 
m p ) for g /A = 2 when the backreaction effect is taken into 
account. The horizontal line corresponds to the Toda-Brumer 
border. 



backreaction effect becomes important around x = 60 in 
this case. 

By implementing the backreaction as a Hartree ap- 
proximation, we find that this typically tends to work to 
suppress the growth of field fluctuations. As illustrated in 
Fig. [7] the fluctuations do not exhibit rapid growth after 
the backreaction begins to work (x > 60). Nevertheless 
we need to caution that linear perturbation theory is no 
longer valid at this stage. For completeness it is required 
to account for the mode-mode coupling (rescattering) be- 
tween the fluctuations @. In fact we found a numerical 
instability for x > 80 in the simulation of Fig. {7\ which 
signals the limitation of the Hartree approximation. It 
is of interest to sec the effect of chaos at this fully non- 
linear stage, but this is a non-trivial problem because of 
the complex nature of reheating. Note that the chaotic 
period ends at some time to complete reheating. It is 
difficult to judge when chaos ends in our system, since 
the mechanism for the decay of <f> and x after preheating 
is not completely known. 

When g 2 /X = 0(1) one can find out the existence of 
chaos during a short period before the backreaction be- 
gins to work. As we already mentioned, the criterion for 
chaos is given by Eq. I|33|) . whereas the criterion for the 
backreaction corresponds to Eq. (|38() . The initial value 
of the quasi-homogeneous field x S e ^ s smaller for larger 
g 2 /X, as illustrated in Fig. El This means that the condi- 
tion (|38|l tends to be satisfied prior to the time at which 
the quasi-homogeneous field x grows to satisfy the Toda- 
Brumer test l|33|) . In Fig. [H] we plot the evolution of the 
system for g 2 /X = 5000 with the backreaction effect of 
created particles. In this case the backreaction becomes 
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FIG. 8: As in Fig. □ for g 2 /X = 5000. 



important around x = 50 before the quasi-homogeneous 
field x increases sufficiently to satisfy the Toda-Brumer 
test. Therefore, when g 2 /X ^ 1, we do not find a signa- 
ture of chaos before the perturbations reach a nonlinear 
regime. 

For the coupling g that belongs to the resonance bands 
n(2n — 1) < g 2 /X < n(2n + 1), we find that chaos occurs 
for 

g 2 /X < O(10) , (39) 

before the backreaction begins to work. We note that 
this result is obtained by using the value J (5x 2 (tf))k<k f 
derived in Sec. II as the initial condition of x at the be- 
ginning of preheating. When g 2 /X > O(10) the field 
fluctuation (Sx 2 ) satisfies the condition (|38[) prior to the 
time at which the necessary condition for chaos is ful- 
filled. Provided that g 2 /X 3> 1, the standard Floquet 
theory of parametric resonance |3lLl32 | is valid at the lin- 
ear level, since the effect of the quasi-homogeneous field 
X is not important relative to its perturbations on sub- 
Hubble scales. 



V. QUADRATIC POTENTIAL 

For the quadratic potential the system can not be re- 
duced to the analysis in Minkowski spacctime by intro- 
ducing conformal quantities. Therefore the analysis in 
the quadratic potential is more complicated than in the 
case of the quartic potential in an expanding background. 

We start our analysis by studying the two-field dynam- 
ics in a frictionless background. In this case the fields 
oscillate coherently without an adiabatic damping due 
to cosmic expansion. The system has the field equations 



corresponding to H — in Eqs. J5J) and (J3J) together with 
the constraint 

E=\4> 2 + \x 2 + V{^x), (40) 

where E is conserved. Unlike the case of an expand- 
ing background, the field x is enhanced only when the 
system is inside of resonance bands from the beginning 
of preheating. This comes from the fact that the field 
X does not shift to other stability /instability bands in 
the absence of cosmic expansion We wish to study 
the existence of chaos for the coupling g that leads to 
parametric excitation of x m an expanding background 
(g > 3.0 x 10 -4 ). First we carry out the analysis in a con- 
served Hamiltonian system given above and then proceed 
to the case in which the expansion of universe is taken 
into account. 

As illustrated in Fig.^ the quasi-homogeneous x field 
at the end of inflation is estimated to be Xf ~ 10 -9 - 
10~ 8 ?Tip for the coupling g > 3.0 x 10 -4 , which is smaller 
than that in the quartic potential with g 2 /X = 0(1). 
Hence the field x f° r the quadratic potential is more 
strongly suppressed during inflation for the values of g 
relevant to efficient preheating. 

In Fig. |5| we plot the evolution of the background field 
X together with y/ {Sx 2 ) and yj (54> 2 ) for the coupling 
g = 3.0 x 10~ 4 . Note that we implement the backre- 
action of sub-Hubblc field fluctuations as a Hartree ap- 
proximation. In this case the Toda-Brumer test gives 
the condition x/ m p > 1-9 x 10~ 3 by Eq. For the 

quadratic potential the backreaction begins to work for 
(^X 2 ) ~ m2 /g 2 , which is basically a similar condition to 
the Toda-Brumer test for the background field x- As 
shown in Fig. [5] the quasi-homogeneous field x does not 
satisfy the Toda-Brumer test, since the backreaction be- 
comes important before x grows sufficiently. 

The quasi-homogeneous x field at the end of inflation 
gets smaller for larger g, see Eq. (|12fl . FigurelTUlshows the 
evolution of the system for g = 3.0 x 1 in which case 
the initial variance of the field x is suppressed relative 
to the case g = 3.0 x 10~ 4 . Although the condition for 
chaos using the Toda-Brumer test l|32[) gets milder for 
larger g, this property is compensated by the suppres- 
sion of the quasi-homogeneous field x at the beginning 
of preheating. Therefore it is difficult to satisfy the nec- 
essary condition for chaos before the backreaction begins 
to work. We carried out numerical simulations for other 
values of g and found that the signature of chaos is not 
seen in the frictionless system as long as the backreaction 
effect is taken into account. 

If we implement the effect of cosmic expansion, the 
energy of the system given by Eq. H4U|) decreases. In fact 
the time-derivative of E is given by 

dE 

— = -3H(cb 2 + x 2 )^-3HE, (41) 

where we used the approximation E ~ <f> 2 + x 2 ■ Then the 
energy lost during one oscillation of the inflaton (At ~ 
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FIG. 9: Evolution of x, (8 X 2 ) 1/2 and {5cj> 2 ) 1/2 (normalized by 
nip) for g = 3.0 x 10~ 4 without the friction due to cosmic 
expansion. The horizontal line represents the border of the 
Toda-Brumer test. We implement the backreaction effect of 
created particles as a Hartree approximation. 



condition, \AE/E\ <C 1, is satisfied. 

In addition to the energy loss, we need to caution 
that the structure of resonance changes in the presence 
of cosmic expansion. The field x passes many instabil- 
ity/stability bands, which is generally called stochastic 
resonance |lOl | . Parametric resonance ends when the res- 
onance parameter, q — g 2 4> 2 /(4m 2 ), drops down to less 
than of order unity, whose property is different from the 
analysis in Minkowski spacetime. 

In spite of above complexities, it is possible to check 
whether the signature of chaos is seen or not in a linear 
perturbation regime. Note that the Toda-Brumer condi- 
tion is valid in an expanding background. We run our 
numerical code including the backreaction of sub- Hubble 
field fluctuations for the coupling g relevant to efficient 
preheating (<? > 3.0 x 10~ 4 ) and find that the background 
X stops growing by the backreaction effect before the 
Toda-Brumer condition is satisfied. This property is sim- 
ilar to the case discussed in the frictionless background. 
Therefore chaos does not appear for the quadratic po- 
tential at least in a regime before the backreaction sets 
in. 



VI. CONCLUSIONS 
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FIG. 10: As in Fig.^with g = 3.0 x 10~ 3 . 



\/m) is estimated as 



(42) 



Therefore the energy loss is large at the beginning of 
preheating (\AE/E\ > 0.1), but it becomes smaller and 
smaller with time. The analysis in the frictionless sys- 
tem can be used in an expanding background when the 



In this paper we discussed chaotic dynamics in 
two-field preheating with monomial inflaton potentials 
V(<p) — Vo(j) n . A scalar field x coupled to inflaton with an 
interaction (l/2)g 2 (j) 2 x 2 is amplified by parametric reso- 
nance when the inflaton oscillates coherently. As long 
as the background field x IS n ot strongly damped in an 
inflationary epoch, x can grow to the same order as 0, 
thereby giving rise to a possibility of chaos during pre- 
heating. 

First we estimated the amplitude of the quasi- 
homogeneous x field at the beginning of reheating by 
considering the variance of the x fluctuations on scales 
larger than the Hubble radius at the end of inflation. 

The variance J (Sxt(tf))k<k f is plotted as a function of 
the coupling g in Figs. Hand El For the quartic inflaton 
potential with g 2 /X = 0(1) there exist large-scale pertur- 
bation modes (k < k c ) which are not strongly suppressed 
during inflation. Then this gives large contribution to 
the total variance of x relative to the case g 2 /X ^> 1 at 
the end of inflation. 

Typically the background field x 1S assumed to be neg- 
ligible in standard analysis of particle creations in pre- 
heating, but its presence leads to a mixing between two 
fields. This gives rise to a chaotic instability in addi- 
tion to the enhancement of perturbations by parametric 
resonance. We note that standard Floquet theory using 
Mathieu or Lame equation ceases to be valid when the 
chaos is present. This new channel of instability can alter 
the maximum size of field fluctuations if two dynamical 
fields do not decay for a long time. 

In order to study whether chaos really appears or not, 
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we solved the background equations J3J|, and to- 
gether with perturbed equations l|23|) and (|24|) . For a 
quartic potential (n = 4 and Vb = A/4) parametric res- 
onance can occur for the coupling g 2 /X of order unity. 
In this case the quasi-homogeneous field x satisfies the 
Toda-Brumer test t[A'A\ before the backreaction effect of 
created particles becomes important. Since this is only 
the necessary condition for chaos, we also evaluated Lya- 
punov exponents for g 2 /A = 0(1) in order to confirm the 
presence of chaos. We find that the maximal Lyapunov 
exponent begins to increase toward a positive constant 
value after the Toda-Brumer condition is satisfied, which 
shows the existence of chaos. Our analysis using a Frac- 
tal map also implies that the system exhibits chaotic be- 
havior when g 2 / X is not too much larger than unity. For 
larger g 2 /A, the backreaction of field fluctuations on sub- 
Hubble scales works earlier than the time at which the 
Toda-Brumer test is satisfied for the background field x- 
We find signatures of chaos for g 2 /A < O(10) at the lin- 
ear regime before the backreaction begins to work. For 
g 2 /X ^> 1 the quasi-homogeneous field x gets smaller at 
the beginning of preheating, in which case perturbations 
enter a nonlinear region before chaos can be seen. 

The system with a quadratic (n = 2) potential can 
not be effectively reduced to the analysis in Minkowski 
spacetime unlike the quartic potential. We first analyze 
the preheating dynamics in a frictionless background and 
then proceed to the case in which the expansion of uni- 
verse is taken into account. In this model the background 
field x does not grow sufficiently to satisfy the Toda- 
Brumer test in the presence of the backreaction effect of 
created particles. We find that this result holds both in 
Minkowski and expanding backgrounds. Therefore chaos 
can not be observed at least at the linear stage of pre- 
heating before the backreaction sets in. 



When chaos is present, this means that the second field 
X gives a non-negligible contribution for the background 
dynamics. This generally gives a strong correlation be- 
tween adiabatic and isocurvature metric perturbations 
|38l |. which can lead to the amplification of curvature 
perturbations for the self-coupling inflation model |3flj . 
This reduces the tensor to scalar ratio r, which is favor- 
able from the observational point of view. Of course the 
ratio r is not sufficient to judge whether the model is 
rescued or not, since large contribution of isocurvature 
perturbations modifies the CMB power spectrum. It is 
still interesting that the appearance of chaos has a pos- 
sibility to reduce the ratio r. 

In this work we did not study the nonlinear dynamics 
of the system during which the mode-mode coupling be- 
tween perturbations plays an important role. We expect 
that the presence of chaos persists even in such a non- 
linear stage provided that two interacting fields <fi and 
X are dynamically important. The chaos would finally 
disappear after the energy densities of scalar fields are 
converted to that of radiation. It is of interest to ex- 
tend our analysis to such a regime including Born decays 
of scalar fields for a complete understanding of chaos in 
reheating. 
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